On Modeling Perfectly Conducting Sharp Corners 
With Magnetically Inert Dielectrics Of Extreme 

Complex Permittivities 

Constantinos A. Valagiannopoulos and Ari Sihvola 
Department of Radio Science and Engineering, 
School of Electrical Engineering, Aalto University, Finland, 

PO Box 13000, FIN-00076 AALTO. 
Email: {konstantinos. valagiannopoulos, ari. sihvola} @aalto.fi 



Abstract — The idea of replacing an edgy perfectly conducting 
boundary by the corresponding interface filled with a dielectric 
material of extreme complex permittivities, is examined in the 
present work. A semi-analytical solution to the corresponding 
boundary value problems is obtained and the merit of the 
modeling has been checked. Certain conclusions for the effect of 
the constituent material parameters and the geometric features 
of the configuration on the model effectiveness, are drawn and 
discussed. 

Index Terms: Boundary conditions, extreme-parameter ma- 
terials, magnetically inert dielectrics, PEC boundaries, sharp 
corners. 

I. Introduction 

Fictitious, non-penetrable boundaries such as perfect elec- 
trically conducting (PEC) surfaces are frequently used in 
electromagnetic problems to simplify the solution process. 
Impedance boundary conditions have been also extensively 
utilized for the same reason in electromagnetic modeling |Q]; 
in particular, an interesting historical review focusing on the 
philosophy of such a concept is given in Q. The condition 
of a perfect magnetic conductor (PMC) is a very widely-used 
concept in electromagnetic modeling, and several designs have 
been presented in the literature to fabricate artificial magnetic 
conductor surfaces OJ, fl4). Moreove, PEC and its dual, PMC, 
can also be seen as limiting cases of a more general ideal 
boundary: the perfect electromagnetic conductor [5], which is 
employed for treating bianisotropic scattering through certain 
transformations. Still a further class of ideal boundaries is 
determined by forcing conditions on the normal components 
of the electric and magnetic flux densities (instead of the 
tangential components of the fields) which leads to the so- 
called DB boundary and its generalizations [6|. 

In the theoretical and computational analysis, when mod- 
eling antennas and other electromagnetic structures, these 
ideal boundary conditions help in confining the computational 
domain into the region of interest and hence reducing the 
required amount of computations. However, once the design 
has potential to be fabricated, these mathematical conceptions 
should be materialized properly in actual structures. The first 
step towards realization of such boundaries involves replacing 
the boundary by an interface with penetrable substances which 



usually needs to be of extreme permittivity and permeability 
as shown for the PEC case in [7j. However, a full realization 
needs to proceed still further: first to find out what type 
of material microstructure reproduces the macroscopically 
extreme parameter values, and once the recipe has been found, 
finally to fabricate it using available materials, like metals and 
dielecttic media. 



Problems of scattering and diffraction by infinite metallic 
wedges have been thoroughly studied due to their canonical 
importance and suitability to analytical modeling. They are 
widely used to model certain configurations as such parts of 
corner reflectors, open-ended waveguides, transmission lines 
and other devices with polygonal shape existing in a laboratory 
of electromagnetics. A rigorous solution to the problem of 
two identical perfectly conducting parallel wedges loaded with 
a dielectric cylinder, is proposed in [8J. In addition, wave 
diffraction by a wedge with anisotropic impedance boundary 
conditions has been explicitly solved through a probabilistic 
random walk method J9). Finally, a sharp PEC corner with 
a cylindrical protecting cap is investigated in ifTUll . where 
possible applications are also presented. 

In this work, we combine the two aforementioned topics 
(boundary modeling, sharp corners) to examine the possibility 
of materializing a PEC wedge with use of lossy dielectric 
media. We consider a two-dimensional configuration consti- 
tuted by two thin consecutive slabs (hands) excited by a plane 
wave. Due to the finite size of the construction, we perform a 
discretization of the hands into a large number of tiny cylinders 
and a semi-analytic solution is obtained. The actual current of 
the metallic case is compared to the polarization current of the 
dielecttic case and similar comparisons have been made for the 
scattering field. Furthermore, the difference in the developed 
field close to the edge, is represented with respect to the 
real and imaginary permittivity of the dielectric, with several 
angular extents and hand lengths, for various slab thicknesses 
and excitation angles. In this sense, certain conclusions on the 
effectiveness of the PEC wedge modeling via magnetically 
inert materials are stated and explained. 
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Fig. 1. The configuration of the comer structure. The two thin consecutive 
slabs forming a wedge can be either metallic or constructed from magnetically 
inert dielectric material. 



II. Mathematical Methodology 
A. Primary Definitions 

In Fig. [T[ the physical configuration of the considered 
construction is shown, and also the interchangeably used co- 
ordinate systems (Cartesian (x,y,z) and cylindrical (p,(f>,z)) 
are defined. Two thin slabs, with different lengths Ly,L?,, are 
forming a corner of angular extent 28, placed symmetrically 
with respect to x axis. The thickness of the strips is not 
infinitesimal, yet small, equaling to 2a << Li+L 2 . The whole 
two-dimensional structure is illuminated by a plane wave, 
of unitary magnitude, whose propagation direction makes an 
angle £ with the positive x semi-axis. The described corner 
could be either metallic or from dielectric with intrinsic pa- 
rameters (e r eo, po)- The scope of the present work is to study 
qualitatively and quantitatively the conditions under which the 
dielectric wedge, close to the sharp corner, responds as similar 
to the metallic analogous as possible. The time dependence 
is of harmonic type e +jwt , where the operating circular 
frequency is denoted by u and the free-space wavenumber 
by k = w^/eoMo = 2tt/A . 

Due to the nature of the configuration and the excitation 
given by: E mc (p,0) = xE inc (p,<p) = ze -^°P«>^-£), the 
electric field in any region will be z-polarized (E(p, <fi) = 
zE(p, (/))). The Green's function, equaling a line source re- 
sponse of electric current j/(upo) positioned along the axis 
(p, (j>) = (P, $), is given by: 
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where J n ,Hn ' are the Bessel and second-type Hankel func- 
tions of rt-th order respectively. To render our boundary value 
problem susceptible to semi-analytical treatment, we regard 
the two hands of the strip wedge separated in a large number 
of consecutive thin cylinders with radius a. Two numberings 



of these pins (u q = 1, ■ 



with are 



to. The positions of the centers of the pins (in the cylindrical 
coordinate system) constituting the upper (first) and the lower 
(second) hand of the wedge, are given by: 

(pi(ui),0i(ui)) = (2aiti,0), ui = 1, ••• ,Ui, (2a) 

(P2(u 2 ), <t>2(u 2 )) = (2au 2 , -0), u 2 = !,-■■ ,U 2 . (2b) 

Finally, the circular boundaries of the pins and the correspond- 
ing areal disks are denoted by c q (u q ) and s q (u q ) respectively 
(q = \,2 , u q = !,-■■ ,U q ). 

B. Metallic Corner 

Let us first consider the metallic wedge. It is well-known 
that a rudimentary analytic tool for solving wave diffraction by 
PEC scatterers, is the radiation integral ifTTl . According to this 
important formula, the scattering field is given by (|3]l which 
is appeared on the top of the next page. The notation k q {u q ) 
corresponds to the z-directed unknown line current flowing 
in the u g -th cylinder. The length variable C is used for the 
line integrations around the circular boundaries of the pins 
with respect to the arguments (Pc,&c)- Due to the small 
thickness 2a of the slabs, we are going to impose the boundary 
conditions for vanishing field around the metallic rods, only 
on (Ui + U 2 ) specific points: the centers of the circular 
bounds: E inc (p q (v q ),(j) q (v q )) + E scat (p q (v q ),<p q (v q )) = 
(u q = 1, • • • , U q , q — 1, 2). That yields to: 
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The vectors of the surface current k g , are containing the 
unknown quantities k q (u q ), while the elements of the matrices 
and the constant vectors are given explicitly in the expressions 
<(3j of the next page. The indexes take the first two positive 
integer values: p,q — 1,2. Once the unknown parameters k q 
are determined, the scattered field into the vacuum background 
can be directly evaluated through (|3]l and the actual induced 
current K act (d) on the metallic wedge can be found as a 
function of the tangential distance on the strip d. 

C. Dielectric Corner 

As far as the dielectric case is concerned, we employ another 
pivotal formula known as scattering integral [12] which, in 
the present case, can be written as in d6j, shown in the next 
page. The notation e q (u q ) has been reserved for the unknown 
electric field on the u g -th circular disk. The surface variable 
S is used for the area integrations on the cross sections of 
the pins with respect to the arguments (Ps,§s)- By testing 
(|6]l on the centers of the consecutive dielectric cylinders, the 
following linear system is obtained: 
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used in proportion to which hand the cylinders are referred 



where I is the (Ui + U 2 ) x (Ui + U 2 ) identity matrix. The 
related explicit forms for the matrices of the dielectric problem 
are given in ([8]) which are also explicitly presented in the next 
page. 

Similarly to the metallic case, the scattered field by the 
dielectric corner is found via ([6]), based on the fields e,j 
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determinable from the system. The axial polarization current 
imitating the sources, producing the scattering field, is defined 
by El: 



K po i{d) = — - — (e r 



l)E(d), 



(9) 



where d is the distance on the strip and E(d) the total electric 
field along the slab hands. 

III. Numerical Results 

A. Parameter Ranges 

Prior to proceeding to the numerical simulation and com- 
menting on the produced graphs, we present the value ranges 
of the input parameters are as follows. Firstly, the operating 
wavelength Ao is kept constant and equal to 1 m (ko — 2ir 
rad/m) throughout the numerical simulations, while the lengths 
of the wedge hands are mainly chosen within the interval: 
< Li,L2 < 2Ao- In this sense, various cases of different 
electrical sizes and symmetries are examined, without caring 
about realistic operating frequencies. The angular extent of 
the wedge covers all possible values (0° < 29 < 180°), while 
the plane wave usually travels along the positive x semi-axis 
(£ = 0°) to excite the sharp edge. When it comes to the 
complex relative permittivity e r = 5R[e, r ] + j!J[e r ], its real 
part can possess extremely large positive values, while the 
opposite happens for the imaginary part due to the selected 
time dependence (SR[e r ] >> l,3[e r ] << —1). The thickness 
of the slabs should be much smaller compared to the hand 
lengths, so that our assumption for thin-cylinder discretization 
is justified, namely a < 0.02Ao. The quantities that are 
of interest include the surface currents K act ,K po i and the 



scattering field E scat . The average relative error in the field, 
around a small cylinder of radius po in the vicinity of the 
edge is expressing how well the dielectric surface imitates the 
metallic one and defined as follows: 



\Em(po,4>) 



DE(po) = 
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(10) 



The symbols Em, Ed are used for the electric fields in the 
metallic and dielectric case respectively, while the integration 
is carried out numerically. 



B. Comparisons and Validations 

One of the scopes of the present report is to examine how 
well a sharp PEC boundary can be materialized with use of 
a magnetically inert dielectric. Given the fact that we are 
mainly interested in the variation of the waves close to the 
edge, it is meaningful to compare the response of our finite 
metallic construction with the case of a PEC infinite wedge. 
Only if the field behavior in the vicinity of the sharp corner 
for the finite and the infinite case is similar, it is worth to 
study the materialization of the sharp PEC surface with least 
penetrable substances. It is well-known lfl4ll . that the solution 
to the plane wave scattering by an infinite PEC wedge with 
the same characteristics as the examined one (angular extent 
29, same placement of the coordinate system as in Fig. [T| is 
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Fig. 2. The: (a) real part and (b) imaginary part of the electric field 
as functions of the azimuthal angle for various lengths of the hands. Plot 
parameters: A = 1 m, B = 45°, £ = 0°, a/X = 0.015, pa = 1 m. 
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Fig. 3. The average field difference compared with the infinite metallic edge 
structure as function of: (a) the radial observation distance (8 = 45°) and 
(b) the angular extent of the wedge for various lengths of the wedge hands 
(po = 0.75A ). Plot parameters: A = 1 m, £ = 0°, a/\ = 0.015. 



given by: 
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with /i m = 2 U-6) • Similarly to ( 10 1, one can define a quantity 
that expresses the difference between the response of the 
infinite and the finite PEC structure: 
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Only if AE(p ) is relatively low, one can proceed in in- 
vestigating the conditions under which a mathematical PEC 
boundary can be replaced by a dielectric interface. 

In Figs [2] we depict the real and the imaginary part of the 
electric field as functions of the azimuthal angle <fi 6 (6, 2ir—9) 
for metallic corners of different Li = L%. It is natural that for 
larger hand lengths the curves of Em (po , <t>) gets closer to 
the solid black line given by Ew(po, 4>)- It is worth noting 



that when the size of the segments equals a couple of free 
space wavelengths, both parts of the electric field are almost 



identical to those of ( 11 1. In Fig. 3(a) the difference AE(po) 
is represented as function of the electrical distance po/Xq for 
several lengths L\ = L2- Again, one notices that the larger are 
the hands of a finite wedge, the closer is the exhibited behavior 
to that of an infinite construction. Moreover, the curves are 
upward sloping with respect to the radial observation distance 
which remarks that the field in the vicinity of the sharp 
edge is not affected by the size of the corner. In Fig. |3(b)| 
the difference in the field between the finite and the infinite 
case is shown with respect to the half angular extent 0. It is 
noteworthy that the more oblique is the corner, the smaller is 
the recorded error, a feature which is justified by the more 
powerful response of the edge. As a conclusion, one could 
point out that the field near the edge (po << Li,^) of an 
infinite PEC wedge is almost equal to that of a finite wedge 
provided that its hands are electrically large (Li,Zv2 > ^o). 

To provide an additional validation check for our techniques 
and computations, we simulated the metallic structure and the 
excitation depicted in Fig. [T] using COMSOL Multiphysics 
commercial software lfT31 . The magnitude of the developed 
(total) electric field at a specific point (p, <f) — (po, cj>o) is rep- 
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Fig. 4. The magnitude of the total electric field at a specific point as function 
of the electrical length of the two equal hands for: (a) 8 = 45° and (b) 
6 = 60°. Plot parameters: A = 1 m, £ = 0°, a/X = 0.015, <j> = 180°, 
PoMo = 0.3. 
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Fig. 5. The actual and the polarization current on the metallic and the 
dielectric obstacle respectively, as functions of the tangential distance on the 
strip for: (a) £ = 0° and (b) £ = 180°. Plot parameters: Ao = 1 m, Li/Xo = 
1, L 2 /X = 1, «r = 800, 6 = 45°, a/X = 0.015. 



resented as function of the electrical lengths Li/Xq = L2/X0 
of the hands which are of equal size. The simulation results 
are compared against those derived from expressions ([3]) and 
(fTT)i for the finite and the infinite wedge respectively. In Fig. 



(Li = L2 = Ao). Regardless of the incidence angle (£ = 0° 



4(a) the angular extent is taken equal to 9 — 45°; it is 
clear that the coincidence between the two curves (blue for 
semi-analytic approach, red for COMSOL numerical method) 
describing the finite wedge is quite satisfying. Naturally, the 
field solution for the infinite wedge (black curve) is constant 
and not dependent on the length of the hands. Mind that 
when the lengths L\ = L 2 are increasing, the difference 
between the field of the finite and the infinite structure gets 



diminished. In Fig. 4(b) where the angle of the wedge is larger 
(6 = 60°), the agreement between the COMSOL simulation 
and the analytical evaluation is also remarkable. 

C. Currents and Fields 

In Figs [5] we show the variation of the actual current K act 
induced on the metallic wedge and the respective polarization 
current K po i, when e r = 800, as function of the normalized 
tangential distance on the strip d/X whose negative values 
correspond to the lower hand of the wedge. The symmetry of 
the graphs is dictated by the related symmetry of the structure 



in Fig. 5(a) and £ = 180° in Fig. 5(a) 1, there is a remarkable 
coincidence between the two quantities which shows that the 
limit of the polarization current for 5J[e r ] — » +00, equals 
the actual current on the PEC surface. When the structure is 
illuminated from the left side (£ = 0°), a maximum current is 
recorded on the corner (d/Xo — 0), while in the case of right- 
sided incidence (£ = 180°), the quantity vanishes at the same 
point. Note also the substantial charge concentration along the 
ending edges of the two strips (d/Xo — ±1), is considerable 
for both excitation waves. 

In Figs [6] the real and imaginary parts of the scattering 
electric component, are represented around a circular bound 
including the whole structure, as functions of the azimuthal 
angle <fi, for the metallic strip and various dielectric slabs of 
different permittivities. Naturally, the best fitting to the PEC 
case (solid black line) is achieved when both 3?[e r ],3[e r ] are 
extremely high (blue dots). Nevertheless, when only one part 
of the complex permittivity is substantial, the coincidence 
of the curves is more satisfying for the opposite part of 
the scattering field. In particular, the red dashed line in Fig. 



6(a) differs much from the black solid one compared to the 



variation of the green dashed one. The opposite behavior is 
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Fig. 6. The distribution of the: (a) real part and (b) imaginary part of the 
scattered electric field as functions of the azimuthal angle for various relative 
complex permittivities of the dielectric. Plot parameters: Aq = 1 m, L\ = 1 
m, L 2 Ao = 0.5, 9 = 30°, £ = 0°, a/X = 0.015, p /\ = 1.2. 



exhibited in Fig. |6(b)| 

D. Average Errors 

In Figs [7] the average relative error DE(p ) is depicted (in 
logarithmic plot) for several lengths of the strip hands, with 
respect to the real and the imaginary part of the permittivity 
e r . The error gets much more diminished for larger real parts 
than for more substantial imaginary parts. In other words, it 
is preferable to model a PEC wedge with a strong dielectric 
instead of using an equally strong conductor. Note also the 



resonance oscillations of Fig. 7(a) recorded for somehow 
small 5Ji [e r ], in contrast to the smoothly decaying curves of Fig. 
|7(b)| Mind that the three curves in Fig. |7(b)| (and partially in 
Fig. |7(a)) , are almost identical each other; this indicates that 
the lengths of the wedges's hands do not play a significant 
role in the behavior of the device which is dominated by the 
presence of the sharp corner. However, the most fascinating 
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Fig. 7. The average relative error of the electric field around the edge 
as function of: (a) the real part of the relative complex permittivity of 
the dielectric and (b) the absolute imaginary part of the relative complex 
permittivity of the dielectric, for various lengths of the lower hand. Plot 
parameters: A = 1 m, Li/A = 1, 9 = 30°, § = 0°, a/A = 0.015, 
po/Ao = 0.1. 



a critical point of 5R[e r ], above which the electrical thickness of 
the strip gets stronger, namely the transmission mechanism is 
activated and gradually prevails over the competing reflection 
mechanism. In other words, the reflection is closer to 100% 
(characteristic of the metallic surfaces) when the slab is both 
electrically dense (large 5R[e r ]) and electrically thin (small 
2koay/^R[e r ]), otherwise it behaves more as a layer than a strip 
and therefore the transmission through itself is not negligible 
("large-a limit"). 

In Figs [8] the average relative error DE(po) is shown as 
function of both parts of e r , for various angular extents of the 
considered wedge. The results are similar to those of Figs [7] as 
the aforementioned minimum ("large-a limit"), where the relay 
from reflection to transmission mechanism happens, is also 
noticed in Fig 



8(a) Nonetheless, the curves in Fig. 8(b) 



more dispersive compared to the the respective of Fig. 7(b) 



are 
a 



feature of these results concerns the minimum of Fig. 7(a) 
observed close to 5ft [e r ] = 450. It seems rather peculiar, given 
the fact that the larger permittivity is used, the smaller the 
error should be measured. Such a seeming contradiction is 
justified by the finite thickness of the wedge; in fact, there is 



feature that remarks the effect of the corner angle (much more 
than this of the hand lengths) on the response of the regarded 
structure. Again the average errors for increasing real part in 



Fig. 7(a) are lower than those of Fig. 7(b) a feature that is 
attributed to thermal losses induced by the nonzero material 
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Fig. 8. The average relative error of the electric field around the edge 
as function of: (a) the real part of the relative complex permittivity of 
the dielectric and (b) the absolute imaginary part of the relative complex 
permittivity of the dielectric, for various angular extent of the corner. Plot 
parameters: Arj = 1 m, Li/Xq = 1, L2/X0 = 1, £ = 0°, a/Xo = 0.015, 
po/Xo = 0.1. 
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Fig. 9. The average relative error of the electric field around the edge as 
function of: (a) the real part of the relative complex permittivity of the dielec- 
tric and (b) the absolute imaginary part of the relative complex permittivity 
of the dielectric, for various thicknesses of the strip. Plot parameters: Ao = 1 
m, Li/A = 1, L 2 /A = 1,9 = 45°, £ = 0°, p /A = 0.1. 



conductivity. 

In Figs [9] the error quantity DE(po) is represented with 
respect to real and imaginary parts of e r for several thick- 
nesses a/Ao of the strips. The described transition from one 
mechanism to the other for varying 5ft[e r ], is clearly shown 
in Fig. |9(a)| where the curves reach their minimum levels at 
larger 3?[e r ] for thinner slabs. In Fig. 9(b)| one observes that 
when the thickness is small, better modeling is achieved with 
larger |3[e r ]|, when the structure behaves as a thin conducting 
strip. On the contrary, more accurate results are recorded when 
considering thicker slabs with small-|3[e r ]| dielectrics when 
the structure reacts chiefly as a rectangular scatterer. 



IV. Conclusions 

An attempt to imitate the reaction of a sharp metallic corner 
to an incident wave with use of magnetically inert materials, 
has been performed in this study. A semi-analytic technique 
is developed to compare the scattering field in both cases. 
Several graphs of the field difference between the metallic and 
the dielectric case with respect to the complex permittivity of 



the constituent material and the geometrical characteristics of 
the structure, are presented and interpreted. 

An interesting expansion of the described work would be 
to use magnetically active and anisotropic materials to model 
sharp metallic boundaries. Additionally, different excitations 
could be considered in order to test the robustness of the 
proposed PEC modeling with respect to different types of 
electromagnetic interaction. Finally, similar approaches would 
be useful in modeling other types of boundaries such as perfect 
electromagnetic wedges or sharp impedance surfaces. 
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